DITHERED TIME DELAYS AND CHAOS SYNCHRONIZATION IN LASERS WITH 

MULTIPLE FEEDBACKS 
E.M.Shahverdiev ^'^ and K. A. Shore ^ 
^School of Electronic Engineering, Bangor University, Dean St., Bangor, LL57 lUT, Wales, UK 
^Institute of Physics, H.Javid Avenue, 33, Baku, AZ1143, Azerbaijan 

ABSTRACT 

By studying laser systems with multiple time delays, we demonstrate that the signatures of time 
delays in the autocorrelation coefficient and the mutual information of the laser output can be 
erased for systems with variable time delays. This property makes such laser systems highly suit- 
able for secure chaos-based communication systems. We also present the first report on chaos 
synchronization in both unidirectionally and bidirectionally coupled variable multiple time delay 
laser diodes with electro-optical feedbacks. 

PACS number(s):42.55.Px, 42.65. Sf, 05.45.Xt, 42.60.Mi,05.45.Gg, 05.45.Vx 

Introduction. -In recent years due to its fundamental and applied interests chaos synchronization 
has been the subject of extensive studies [1-2]. From the application point of view, chaos based 
communication systems can offer improved privacy and security in data transmission, especially 
after the recent field demonstration using a metropolitan fibre network [3] . In optical chaos based 
communications, the chaotic waveform is generated by using semiconductor lasers with all-optical 
or electro-optical feedback(s). Mathematically such time delayed systems are described by delay 
differential equations (DDE) [4]. Recent studies have revealed that DDE-chaos based communi- 
cation systems [1-2] could be more secure than one based on chaos in ODE (ordinary differential 
equations) [5], as DDEs usually have a very high dimensionality and many positive Lyapunov ex- 
ponents. However, the authors of [6], were successful in extracting messages masked by a chaotic 
signal of a time delay system. In [6] the authors used the average fitting error as a function of 
the embedding delay, which had a pronounced minimum at the time delay. This allowed them 
to correctly identify the time delay of the system and to extract the message successfully using 
a simple local reconstruction of the time delay system. The delay time can also be exposed by 
several other measures, e.g. autocorrelation coefficient [7] and mutual information [8]. 

Quite recently it was suggested that multiple time delay systems can offer a higher complexity 
of dynamics than is achievable in more conventional single delay time systems [9]. The role of 
additional time delays in achieving a homogeneous steady state in coupled chaotic maps was inves- 
tigated in recent work [10] with the stabilizing role of the additional time delays being emphasized. 
The stabilizing role of additional time delays in external cavity laser diodes with multiple feedbacks 
is highlighted in [1 1] . The high complexity of such multiple time delay systems can provide a new 
architecture for enhancing message security in chaos based encryption systems. However using 
the same approach as in [6], it was shown that for the multiple feedback systems, time delays can 
also be successfully recovered and reconstruction of the dynamics of the systems was still possible 
[12]. A viable recent proposal to enhance the security of the chaos-based communication schemes 
is the use of time delayed systems in which the delay time is modulated in time [13]. 
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In this paper we study semiconductor lasers with modulated multiple time delay electrooptical 
feedbacks. First we show that contrary to the case of constant time delays, with variable time 
delay systems the delays cannot be recovered by investigating the autocorrelation coefficient or 
the mutual information of the laser output. Thus variable time delay systems can offer an increase 
in the security of chaos based communication systems. For message decoding in such schemes one 
has to be able to synchronize the transmitter and receiver lasers. In the light of this, in the paper 
we also present the first report on chaos synchronization in both unidirectionally and bidirection- 
ally coupled variable multiple time delayed semiconductor lasers with electro-optical feedbacks. 
Wavelength chaos model.-The laser system considered in this paper is an electrically tunable Dis- 
tributed Bragg Reflector (DBR) laser diode with a feedback loops, whose property is to exhibit 
a nonlincarity in wavelength. This system was proposed in [14] as a chaotic wavelength signal 
generator for chaos based secure communication. 

The output wavelength of the chaotic oscillator is described by the following dynamical equa- 
tion: = -n't) + /5Asin^(-^A(t - r) - $0), where A is the wavelength deviation from the 
center wavelength Aq;D is the optical path difference of the birefringent plate that provides the 
nonlincarity; $0 is the feedback phase;r-the feedback loop delay time;T is the time response in the 
feedback \oop;f3x is the feedback strength. With x = {7iDX)Aq'^ and m = {t:D^x)A.q'^ we rewrite 
it in the following normalized form:^^ = —ax{t) + msin^(a;^ — $0)) where a is the relaxation 
coefficient and Xt = x{t — r). In this paper we will consider the variable time delay case r = r(t). 
In the following we consider chaos synchronization between both unidirectionally and bidirection- 
ally coupled identical laser diodes with double feedback. In figure 1 a schematic diagram of the 
experimental set-up for synchronization is given. 
The principal elements of the chaotic oscillator are: 

(1) a tunable laser source (DBR) whose wavelength can be tuned continuously, i.e. by a DBR- 
section injection current /; 

(2) a wavelength nonhnear element formed by a birefringent plate (BP) set between two crossed 
polarizers, inducing sin^-nonlincarity; 

(3) a photodetector (PD) providing a linear conversion of the optical power into a photocurrent 
with a conversion factor G, which can be adjusted electronically; 

(4) a delay line (DL) which introduces a time delay much longer than the response time of the 
feedback loop, in order to obtain the chaotic regime; 

(5) a first-order low-pass filter (LPF), which determines the response time of the feedback loop. 
Generally, an experimental realization of the variable time delays can be achieved by changing of 
the distance between the light source and photodetector, e.g. by using trombone structure, as for 
example in [15]. 

Numerical simulations. -First we consider the case of unidirectionally coupled laser systems: The 
master laser is described by the following equation 




(1) 



The dynamics of the slave laser is governed by the following equation: 



= -a2y{t) + 1713 sin^iVr, - b<^o) + sin^y^^ - $o) + sin2(x^3 - $o), (2) 

where ti,2 = toi,o2 + xi{t)Tai,a2 sin{u!i^2t) are the variable feedback loop delay times;T3 = tqs + 
X 1(^)^03 sin{uj3t) is the variable time of flight between master laser x and slave laser y;Toi,o2,o3 are 
the zero- frequency component, rai a2,a3 are the amplitude, ^^'^'^ are the frequency of the modula- 
tions; Xi(t) is the output of laser (1) for constant time delays, i.e. Ti = Toi,T2 = T02; Tni^2 and 7713 4 
are the feedback strengths for the master and slave laser, respectively; K is the coupling strength 
between the lasers. For mutually coupled systems the term Ksm'^{xr3 — $0) should be added to 
the right-hand side of Eq.(l). In the experimental scheme unidirectional coupling can be realized 
by inclusion of an optical isolator (01), figure 1. 

In the case of variable time delays establishing the existence and stability conditions for the 
synchronization is not as straightforward as for the constant time delays. Having in mind that 
for a; = we obtain a case of constant time delays, then as an initial guess one can benefit 
from the existence conditions for the constant time delays case, see,e.g.[16]. It is our conjecture 
that high quality complete synchronization x = y will be obtained if the parameters satisfy the 
conditionsimi = 7713 + 7^, m2 = 7714 if Ti(t) = T^lt) (i.e. the zero frequency components, amplitudes, 
and modulation frequencies are the same:roi = tqs, r^i = Tas, ooi = W3),or 7712 = 7774 + K, nii = 7773 
if r2(t) = 7'3(t)(i.e.ro2 = To3,Ta2 = Ta3,oj2 = i^s)- As evidenced by the numerical simulations below, 
this conjecture is found to be well-based. 

Before studying the synchronization between the laser systems with variable time delays we in- 
vestigate the autocorrelation coefficient and mutual information for the output of the master 
laser for both constant and variable time delays. The autocorrelation coefficient Ca is a measure 
of how well a signal matches a time shifted version of itself and is a special case of the cross- 
correlation coefficient [17] C{At) =< {x{t)- < x >){y{t + At)- < y >) > {< {x{t)- < x >Y >< 
{y{t + At)— < y >Y >)~°"^ for x = y : where x and y arc the outputs of the interacting laser 
systems; the brackets< . > represent the time average; At is a time shift between laser outputs. 
This coefficient indicates the quality of synchronization:C=0 implies no synchronization; C = ±1 
means perfect (inverse) synchronization. The mutual information J(r) between x and Xr is de- 
fined [8] by J(t) = Ex(t),x{t-r)P{^{i),^{i - T))^og2ip{x{t),x{t - T)){p{x{t))p{x{t - t)))-\ where 
p{x{t),x{t — t)) is the point probability and p{x(t)) and p{x{t — r)) are the marginal probability 
densities. The mutual information J measures the information shared by two variables x{t) and 
x(t — r), namely it measures how knowledge of one of these variables reduces the uncertainty 
about the other. 

Constant and variable time delay systems 

Figures 2(a) and 2(b) demonstrate the autocorrelation coefficient and the mutual information for 
the output of laser x for constant time delays, respectively, i.e. for Ui = 002 = 0, «i = «2 = 4, $0 = 
7r/4, 7?7i = 12, 7772 = 15, Tqi = 3, To2 = 5. It is clearly seen that time delays can be easily recovered 
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from both the autocorrelation coefficient and mutual information, as they exhibit extrema at time 
delays or their multiples and combinations. It is noted that the identification of delays in the case 
of two delayed feedbacks was reported in [18]. 

Next let us consider the variable time delays scenario. In investigating the behavior of the 
autocorrelation coefficient and the mutual information we have experimented with different types 
of variable time delays, among them:(a)sinusoidal modulations:ri^2 = toi,o2 + Tai,a2 sin((^i^2^); 
(b)chaotic modulations:ri^2 = Toi,o2 + 'Tai,a2Xi{t);a.nd (c)combined chaotic and sinusoidal modu- 
lations ri,2 = Toi,o2 + xi{t)Tai,a2 sm{u!i^2t) ■ Extcusive numerical simulations have established that 
erasure of the signatures of time delays in the autocorrelation coefficient and the mutual informa- 
tion is best achieved for combined chaotic and sinusoidal modulations of T{t). Figures 3(a) and 
3(b) show the autocorrelation coefficient and the mutual information of the master laser output 
for sinusoidal modulations for Ti{t) = 3 + 0.03 sin(0.006t) and r2(t) = 5 + 0.03 sin(0.006t) with 
other parameters as in figure 2. In this connection we emphasize that the retrieval of the per- 
oidic time delay from experimental time series by use of the mutual information and modified 
filling factor was presented in [19]. The autocorrelation coefficient and the mutual information 
of the laser output for chaotic time delays for ri(t) = 3 -|- 0.03x1 (t) and T2{t) = 5 + 0.03x1 (t) 
are portrayed in figure 4(a) and 4(b). Here Xi(t) is the x laser output for constant time delays 
for parameters as in figure 2. Finally figures 5(a) and 5(b) depict the autocorrelation coefficient 
and the mutual information of the laser output for combined sinusoidal and chaotic time delays, 
Eq.(3) for Ti(i) = 3 + 0.03xi(i) sin(0.006i) and T2{t) = 5 + 0.03xi(i) sin(0.006i) with the rest of 
parameters as for figure 2. 

Here wc also would like to emphasize the following point. As noted in [20], with increasing the 
feedback rate it could become difficult or even impossible to recover the time delay from the time 
series. With this in mind we have conducted numerical simulations with both low and high level 
feedback rates. The numerical simulations have demonstrated that with the right choice of the 
sampling rate of the time series one can still retrieve the time delays even for the case of increased 
feedback rates. We also note that the feature underlined in [20] is not generic. As demonstrated in 
[21] for the semiconductor laser with a single optical feedback low feedback rates are more suitable 
than the higher feedback to make the delay indentification difficult. 

Thus, it is evident that combined chaotic and sinusoidal modulations of time delays is most suc- 
cessful in eliminating the signatures of the time delays in the autocorrelation coefficient and the 
mutual information of the laser output. In other words modulation of the delay times gives rise 
to the loss of their signature in the autocorrelation coefficient and the mutual information, and 
therefore can improve the security of chaos based communication systems. 
Synchronization properties 

As mentioned above, in chaos based communication schemes synchronization between the trans- 
mitter and receiver lasers are vital for message decoding. With this in mind we present here the 
first report of chaos synchronization between variable time-delay lasers. In the literature there are 
a few papers on chaos synchronization between chaotic systems with variable time delays [13,22]; 
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however those papers dealt with synchronization in maps, or simple systems, such as Lorenz, 
Rossler and Mackey-Glass. These paradigm chaos models in nonlinear dynamics are of limited 
practical interest in applications in fast communication schemes. 

Figures 6(a) and 6(b) portray complete chaos synchronization x = y between unidirectionally cou- 
pled lasers, Eqs.(l) and (2) for variable feedback time delays Ti{t) = 3 + 0.03sin(0.006i), T2(t) = 
5 + 0.03sin(0.006i) and variable coupling time delay Ti{t) = 3 + 0.03 sin(0.006t) with parameter 
values as a\ = a2 = 4, mi = 12, m2 = m4 = 15, ma = 0.5, = 11.5,$ = 7r/4. Figures 7(a) 
and 7(b) depicts chaos synchronization between mutually coupled lasers x and y with the same 
variable feedback and coupling time delays as in figure 6 for parameters cti = 0^2 = 4, mi = ma = 
12, m2 = rrii — 15, X = 12.5,$ = 7r/4.(We recall that in the case of mutually coupled lasers 
the the right-hand side of Eq.(l) is augmented by a term i^'sin^(xT-3 ~ ^o)-) The vahics of the 
cross-correlation coefficients for both unidirectionally and bidircctionally coupled systems testify 
to the high quality chaos synchronization, which is vital for information processing in chaos-based 
communication systems. 

In this connection we also study the effect of parameter mismatches on the synchronization 
quality. It is noted that sensitivity of the synchronization to mismatches of the parameters can 
lead to a high level of security due to the difficulty to replicate the receiver laser, i.e. sensitivity 
to parameter mismatches increases the security of encryption. However, the internal parameters 
of the interacting laser diodes unlikely to match exactly even if they are produced from the same 
wafer. Moreover, the operating parameters cannot be perfectly controlled. In other words, in 
practical cases, synchronization must therefore occur also for small parameter mismatches. Most 
importantly, an investigation of the effect of parameter mismatches on synchronization quality 
will enable determination of the most sensitive synchronization parameters. In our numerical sim- 
ulations we allow a 5% mismatch between parameters. 

In figure 8 the dependence of the cross-correlation coefficient C between the unidirectionally cou- 
pled lasers on the ratio ^ of the feedback strength of the transmitter m2 to the feedback strength 
of the receiver m4(<C>), on the ratio ^ of the relaxation coefficient of the transmitter laser ai to the 
relaxation coefficient of the receiver laser q;2(V), on the ratio ^ of the frequency of the feedback 
modulation of the tranmitter uii to the frequency of modulation of the injection term a;3(A), and 
on the ratio ^ of the fixed time delay of the transmitter laser tqi to the fixed time delay of the 
injection term ro3(*) is presented. 

It is emphasized that for these cases the synchronization quality is quite robust for small parameter 
mismatches (1%). These results also demonstrate high sensitivity of the synchronization quality to 
the parameter mismatches of the frequency of feedback modulations and feedback times. We have 
also found the similar trend for the mutually coupled laser systems even with more robustness to 
the parameter mismatches. 

Finally we dwell on the possibility of chaos control via variable time delays. The numerical 
results have demonstrated that there is an optimal frequency of feedback modulations to erase 
the signatures of time delays in the autocorrelation coefficient and the mutual information of the 
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chaotic laser output. We have estabhshed that for frequencies much higher than the optimal one 
the variable feedback(s) can be used to control the chaotic behaviour, i.e. to convert such a be- 
haviour to the fixed state (figure 9). 

Before concluding we also emphasize that the laser model under study is of the Ikeda type and 
can show multi-stability. As it is demonstrated in [23] the relationship between the achievability 
of the synchronization and multi-stability in time delayed systems is a quite interesting subject. 
This relationship deserves more detailed study for the time delayed laser systems which is beyond 
the scope of this paper. 

To summarize, by investigating autocorrelation coefficient and the mutual information we have 
established that in terms of security considerations, variable multiple time delay laser systems 
offer significant advantages for chaos-based communication schemes. We have also reported on 
chaos synchronization in both unidircctionally and bidircctionally coupled variable multiple time 
delay laser diodes. The results of the paper provide the basis for the use of lasers diodes with 
multiple variable time delays in chaos-based secure high-speed communication systems. 
Acknowledgements. -This research was supported by a Marie Curie Action within the 6*'* European 
Community Framework Programme Contract N MIFl-CT-2006-039927 and Contract N MIF2-CT- 
2007-039927-980065. 

Figure captions 

FIG.l. Schematic experimental set-up for wavelength chaos synchronization between the trans- 
mitter and receiver laser diodes with electooptical feedback: DBR LDl and DBR LD2 are the 

Distributed Bragg Reflector transmitter and receiver laser diodes, respectively; BS, beamsplit- 
ter; M, mirror; BP, birefringent plate between crossed polarizers(not shown in the figure); PD, 
photodiode; DL, delay line; LPF, low-pass filter; G, optoelectronic gain; I, DBR-section injection 
current; 01, optical isolator to provide unidirectional coupling between laser diodes. 

FIG. 2. The autocorrelation coefficient Ca (2(a)) and mutual information J (2(b)) of the laser out- 
put for constant time delays, Eq.(l) for ai = a2 = 4, mi = 12, m2 = 15, tqi = 3, T02 = 5, $ = 7r/4. 
Dimensionless units. 

FIG. 3. The autocorrelation coefficient Ca (3(a)) and mutual information J (3(b)) of the laser out- 
put for sinusoidal modulations, Eq.(l) for ri(t) = 3+0.03 sin(0.006i) and T2{t) — 5-1-0.03 sin(0.006i). 
The other parameters are as in figure 2. Dimensionless units. 

FIG.4. The autocorrelation coefficient Ca (4(a)) and mutual information J (4(b)) of the laser 
output for chaotic time delays, Eq.(l) for Ti{t) = 3 -|- 0.03xi{t) and T2{t) = 5 -|- 0.03xi{t). Xi{t) is 
the laser output, Eq.(l) for constant time delays tqi = 3,to2 = 5. The other parameters are as in 
figure 2. Dimensionless units. 
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FIG. 5. The autocorrelation coefficient Ca (5(a)) and mutual information J (5(b)) of the laser 
output for combined sinusoidal and chaotic time delays, Eq.(l) for Ti{t) = 3 + 0.03a;i(t) sin(0.006t) 
and T2{t) ~ 5 + 0.03a;i(t) sin(0.006t). The other parameters are as in figure 2. Dimensionless units. 

FIG. 6. Numerical simulation of unidirectionally coupled variable time delay lasers, Eqs.(l-2) 
for «! = = 4, $ = 7r/4, mi = 12, m2 = 15, ms = 0.5, 777.4 = 15, = 11.5 and Ti{t) = 
3 + 0.03a;i(t)sin(0.006t),T2(t) = 5 + 0.03a;i(t) sin(0.006t), r3(t) = 3 + 0.03x1 (t) sin(0.006t)). Com- 
plete synchronization: 6(a): Time series of the transmitter laser wavelength (x-sohd line) and 
receiver laser wavelength (y-dotted line); 6(b):synchronization error y — x dynamics. G is the 
cross-correlation coefficient between the wavelengths of the transmitter and receiver lasers. Di- 
mensionless units. 



FIG. 7. Numerical simulation of bidirectionally coupled variable time delays lasers, Eqs.(l-2) for 
mi — 7713 — 12,7772 = "7.4 = 15, K — 12.5. The other parameters are as in figure 6. Notice that 
for the mutually coupled variable time delays lasers the coupling term in Eq.(2) should be added 
to the right-hand side of Eq.(l). Complete synchronization: 7(a):Time series of the x laser wave- 
length (solid line) and y laser wavelength (dotted line);7(b):y versus x.C is the cross-correlation 
coefficient between the wavelengths of the transmitter and receiver lasers. Dimensionless units. 

FIG. 8. Numerical simulations of Eqs.(l-2) for unidirectionally coupled lasers. Dependence of 
the cross-correlation coefficient C between the unidirectionally coupled lasers on the ratio — of 
the feedback strength of the transmitter 7772 to the feedback strength of the receiver 7774(0), on 
the ratio — of the relaxation coefficient of the transmitter laser cti to the relaxation coefficient 
of the receiver laser q;2(V), on the ratio ^ of the frequency of the feedback modulation of the 
tranmitter Ui to the frequency of modulation of the injection term Co'3(A), and on the ratio ^ of 
the fixed time delay of the transmitter laser tqi to the fixed time delay of the injection term To3(T*r). 

FIG. 9. Numerical simulation of the variable time delay lasers, Eqs.(l) for 7771 = 12, 7772 = 15, Ti{t) = 
3-1-0. 03^1 (i) sin(lOt) and T2{t) — 5+0.03xi{t) sin(lOi). Chaos control via variable time delays:Time 
series of the x laser wavelength with variable time delays (x(t)-solid line) and with fixed time de- 
lays Toi = 3,To2 = 5 (xi (t)-dotted line). 
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